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SUMMARY 


I"' J^^quency extension of the NASA-Ames unsteady, transonic code 
results created and is evaluated by comparisons with experimental 

Zfiln \ experimental test case is a NACA 64A010 airfoil in pitchinq 

motion at a Mach number of 0.8 over a range of reduced frequencies Com- 

LTRANr^nT’^'’^*^/^'’^ modified code is an improvement of the oViqinal 
ripfjfri Pf'P''iPes closer agreement with experimental lift and moment 
coefficients A discussion of the code modifications, which involve the 

alnoH^n, is’-? to the boundary conditions of the numerical 

aujorithm, is also included. 


I. INTRODUCTION 


°K transonic flow about oscillating airfoils is 

R^earch 0^ bv Ran'ha^s developed at NASA-Ames 

fho f by Rallhaus and Goor.iian in 1976 (ref. 1). LTRAN2 solves 

rPfi nonlinear, sma 1 1 -disturbance equation under the low 

reduced frequency assumption {k z ,a c/U- <0.2, where k is the reduced 
equency bosed on full chord, c). For the purpose of flutter analysis 

calculations indicated the need to perform accurate 

calculations in a frequency range up to k = i.o. With this in mind an 

objective of the present study is to extend the range of reduced frequency 
or improved LTRAN? appl icabi 1 ity. The modifications made to the code in 

the cdlculati3 o7 frequency-time dependent terms in 

the calculation of the pressure coefficient as well as the w.ake and down- 

stream boundary conditions. The low frequency governing equation is re- 

performed similar modifica- 

tnT- neriVaf aioor? r frequency terms to the boundary conditions of 

t e merical algorithm. Their results indicate that under subsonic flow 

anQlP^°7’i agreement with linear theory in amplitudes and phase 

tith V?- "’P^ent coefficients is obtained at higher frequencies 

The mndir. demonstrate the merit of 

e modified^code, calculations done under transonic conditions should be 

compared with experimental results. Hence, this study was undertaken to 
modify the existing code LTRAN2 and to evaluate the effects of these changes 
by comparison with experimental data. 


The test case for comparison (ref. 5) was a NACA fiAAOlO airfoil, pitch- 
inq about quarter chord, at a Mach number of O.R over a ranqe of reduced 
frequencies up to n.6. This case was chosen based on the following cri- 
teria: 

1) the availability of good test data 

2) the presence of a moderate strength shock wave 

3) the absence of strong separation effects in the experiments. 

The numerical procedure including the code modifications used in the 
creation of LTRAN2-HI, a high frequency extension of the original code, is 
discussed in section II. Linear theory comparisons are made in section III 
which indeed confirm the results of other researchers (refs. ?-4) that, in 
general, accuracy is improved at higher reduced frequencies./ Section III 
also presents the results of the nonlinear, experimental comparisons where 
LTRAN2-H1 is shown to be superior. Finally, brief concluding remarks are 
made in section IV, 

The authors wish to thank S. S. Oavis for supplying the linear theory 
results, as well as W. J. McCroskey and G. N. Malcolm for their comments and 
suggestions while reviewing the manuscript. 


II. GOVERNING EQUATION AND BOUNDARY CONDITIONS 


Governing Equation 

An unsteady, transonic smal 1 -disturbance equation in Cartesian coordin- 
ates may be written as 


xt ^xx ^yy (1) 

where ^ ' k^M^^/(S^/’ 

B = k 

c = (I - 

and where is the disturbance velocity potential , Moo is the freestream Mach 

number, y is the ratio of specific heats, and 6 is the airfoil thickness- 
to-chord ratio (ref. fi). The choice of the exponent m is somewhat arbit- 
rary. In the calculations presented here, Spreiter scaling was used (m = 2). 
The parameter k is the reduced frequency. For an airfoil of chord length c, 
traveling with speed ll.v , and executing some unsteady oscillatory motion of 
frequency ui, k z u)c/lKv . Thus the reduced frequency is given in terms of 
radians of oscillatory motion per chord length of airfoil travel. The quan- 



<n equation ( 1 ) have been scaled by c.c/ 6 '/^ ca^/’ li 

; •tranl,r'',Van"d,!^ the’ "fa.'fAar tSoI 

I, transonic, stnal 1 -disturbance equation for steady flows. 

the e*3ea\"o'n’‘'"'"‘'°" <»• 'o'- '<»< freqeencles. la 


2B4' , = C4) + 4 . 

xt XX yy 


( 2 ) 


Lr.^e ‘e" ertVo- '’’d-ee t^raX'^ir" 


k % 6 ^^^ << 1 


( 3 ) 


Boundary Conditions and Code Modifications 
tl™e'^depenTnV\\V 4 ° 7 „ '-TRANJ-H 1 introduce 

T?Sld^"„^tt“ 7 o^Toring°Tqt'’a 7 ic^";Vd°^ Thes?"cPap,?;M'c‘i.’'a"r 7 "^de^! 

discussed below in qre^ate^r detail’. ^ numerical boundary conditions are 
pres^ ^7e7 n c'i ent ^e^Crft ?e7” ’ -^1 Udisterbance 




■2 (v\ + k 4 >j.) 


( 4 ) 


teS'b/Vn 'oTde;TTa^„??rVr:ui”nr‘'''"^[T^irH^l ‘7 ter. was Ct- 

(second ordL ac'ratratVjP' A ''’'wbeVe’T uZ tV?T 

ection) and may be expresselJ: ’ ^ ^ ^ ^ ^ the x dir- 



- 7 - . - ♦j-l \ " 

V AX, .SXj . ; ) 

( - '*'.( • 7 - *.i-i \ 

^ ''7*1, j / 

h"’ - •/ ) ' 



Wake Condition ; 
sheet 


To ensure a continuity of pressure across the vortex 


which now implies 



= n 


[*X * = n (5) 

Where C 3 denote a jump across the wake. The addition of the k-t>f term in 
condition constituted the most significant coding change to 
LTRAN2. Formerly, vorticUy in the wake was assumed to propagate infinitely 
fast downstream (for a characteristic analysis, see reference 1). Circu- 
lation, a function of time alone, was then assumed to be uniform in the x 

dicfoil trailing edge to the downstream boundary, 
with the inclusion of a time dependent term in the pressure coefficient 
calculation, enforcement of the jump conditions for the wake produces a 
circulation quantity now a function of both time and space. Consequently, 
vorticity travels downstream at a finite velocity (in unsealed units, the 

freestream velocity) as shown by the following brief characteristic anal- 
y ^ 1 s * 


circulation r = A^^ , the jump condition from equation 

(5) yields 


k’'t = " 


( 6 ) 


Equation (6) has characteristics 

t = kx + constant 

and : remains constant along these lines. Thus, 1 ’ propagates downstream at 

a vel oc 1 ty 


dx = 1 

dt t 


in scaled units 


Recognizing that 
quant i t i es 


x/c, t 


u)t, k = a)C/Ua> where x, t are unsealed 


dx 


= Ux, 


df. 


I 


Therefore the introduction of the k term in 
accurate description of vorticity propagation 
velocity. 


equation (5) provides a more 
in the wake at the fr^.-estream 


The following Crank -Nicol son algorithm, consistent with the differencing 

froAMo equation, was chosen to implement the wake condition in 

L I K I , 



.. n n 

1 j -1 

AX, 7, 


+ 


n+1 ,, n+1 

AX , , , 


) 




Solving for the value of circulation at the new time level: 


l-.n+l ^ 

J j-1 


n 


’At - 

aTT 

m 


- k Ax 


irlzX 





( 7 ) 


,, Ho»^p stream Boundary Condition ; The downstream boundary 
(formerly = ij as x approaches infinity) becomes 


condi tion 


(}) + k i}>. =0 X » '•’ 

X L 


(8) 


Like the circulation, the disturbance potential 
IS now a function of time and space. 


at the downstream boundary 


rnnHU L '■TRAN2-NLR (ref. ?) downstream boundary 
i™? ^ aeroth order extrapolation (*;, ■- 0) was retained. However, 
implementation of equation (fl) is believed to maintain consistency at th,: 

Specifically. boundary conditions. 


^jniax ^jmax-1 


[ 


At - k Ax . 
j niax. 

At f k Ax.” 


jniax. 


Jniax- 1 
jniax- 1 . 


\^jniax-l 


4>" 

^jtnax 


where jmax is the grid index of the x-location 
equation (7) is also satisfied at j=jniax. 


farthest downstream. 


Hence 


Airfoil Tangency Conditi 
body boundary condition thus 
airfoil tangency condition. 


on: A time dependent term is included in the 

eliminating the low frequency assumption in the 
If .V f(x,t) defines the body surface, then 




(9) 


ensures flow tangency. 


A1 through undocumented, equation (9) (including this time dependent term! 
was used in all published “low frequency" LTRAN2 calculations by Ballhaus 

and 8ooriian. Also, in the calculations presented here, equation (Ql will 
be used by both LTRAN2 and LTRAN2-HI. cMuation will 



m. CALCULATED RESULTS 


Linear Calculations 

Cninpwted results from LTRAN? and LTRAN?-HI are first compared with exact 
T' the capahnit^of the 

thr««tS°,„«r thee"''' ‘"'Ste^y solutions in the linear domain, 

e ac . linear theory results are solutions to equation M) with v - i in 

(?) "als^wi^h T- ^1 to equa'tion 

S^eHcaraldoliZ alnorithm, and ;.) deficiencies in the 

qovernintj equation. 


- •» I ,,s_ I ■ c. I trill, I III i,iit 

resultinq from the npqlect of the term A 4 .^^ in the 

LTRAN?'’wUh'r^inh‘? “f and van der Vooren (ref. ?), 

^ hi qh frequency extension provides better aqreement with linear 

mnmpnt oriqinal LTRAN? in amplitudes and phase anqles of lift and 

vT: r'lce^^ ^ ^ momenfcoeffi 5enfs 

oriqinal I TRAN? provides reasonably accurate results hut only for reduced 

moment"coe%iclent^^ LTRAN^^ '^’^^^^tion of the real component of the 

enc coerncient, LTRAN?-MI provides a more accurate prediction of both 

lifts and moments over the entire ranqe of reduced freruencirs tested 

linear^hIorI;^f^^''’^ comparison of the low frequency LTRAN? with 

Sion is not clelrlC i component of moment), the hiqh frequency exten- 
sion IS not clearly an improvement of the oriqinal code. This fact further 
motivates a comparison with experimental data further 


Nonlinear Experimental Comparisons 

fio previously, comparisons with experimental data for a pitchinq 

iLer rr ref' m’ -«e et , trenso^ic Slch 

nuadv sulutfdo' foJnd’hv ad,iust?np“the ^'oratclirh'’ 

5rr’r!,u!rof''c'''(sTr'iJv)'’- To? 9T'"'’' "ft. mth a 

fomed at a . . -Viioa fherVv ^tch ™ Hm.'"'ru’orprrrr"r^ 

caku at ions were sensitive to the type of smoothing used io defining the 

ta, coorrrrrrvrrVrer:rrs"w;v:^^^^ 

tain some deqree of consistency in model inq the experiment. 

the InsteaTrc^^kuTaHo?^ represent the initial conditions for 

chn.,-? »K- Experimental steady lower surface pressures are 

mt ;^is calculated pressures at the s^e yaTue of 

t. This information should be used as a baseline comparison. 

In the followinq computed unsteady results, both LTRAN? and I tram? m 

IjrrTn" “a'"' ‘‘r for k To maTotain stahilUy 

codes at the lower reduced frequencies, it was necessary to increase 


A 



the number nf time steps per cycle to 96D and 19?0 for k = n.l and k = 0.05, 
respectively. The modifications to LTRAN2 did not adversely affect the sta- 
bility of the code. In fact, LTRAN7-HI, unlike its counterpart, was capable 
of producing a stable solution at k =0.05 with 1440 time steps per cycle 
indicating that the addition of high frequency terms has a stabilizing 
effect on the code. Also, high frequency results may be obtained for essen- 
tially the same "price" as results from the original LTRAN?. 

Figures 4 and 5 display first harmonic comparisons of lift and leading 
edge moment coefficients versus reduced frequency (reference 5 reports first 
harmonic data; computed results were Fourier-analyzed to determine first 
harmonics). Here LTRAN? -HI shows improved agreement with experiment at 
higher frequencies. Note in particular the successful prediction by 
LTRAN2-HI of the critical transition in the leading edge moment from a phase 
lead to a phase lag (fig. 5b). As shown in figures 4a and 5a, the high 
frequency modification, in general , produces an improvement in the calcu- 
lation of real and imaginary components of both lift and moment coefficients 
over the entire reduced frequency range. The greatest improvement is seen 
in the determination of the imaginary components where LTRAN2-HI, unlike the 
original code, captures the experimenta 1 1 y observed trends. For this 
reason, LTRAM2-H! is the recommended version of the code for use in 
transonic r a 1 cul at i ons , 

First harmonic unsteady lower surface pressure comparisons are made in 
figures fi through 1?, spanning the reduced frequency range. Real and 
imaginary components of the pressure coefficient are presented. However, 
amplitude and phase information is included for k = D.l (fig. 7) to display 
the somewhat misleading behavior of toe calculated phase angle after the 
shock (x/c>D.5). In this region the real component of pressure is rel- 
atively small, and therefore a slight error in its computation will 
produce a large change in phase angle ( = tan"' ( '"/Re)). Hence the 
comparison of real and imaginary components with their corresponding 
experimental values is the preferred method of presentation. 

The surface pressure comparisons are found to be inconclusive, however. 
It is difficult to justify the use of I.TRAN2-HI over tne use ot the original 
code based on this information alone, for there are instances where the two 
codes seem equally suitable. But as noted previously in figures 4 and 5, 
LTRAN2-HI provides improvement in lift and moment calculations, especially 
in the prediction of the imaginary components of these loads. Consequently, 
the importance nf examining integrated pressures must not be underestimated. 


IV. CONCLUSIONS 


LTRAN2-HI, a high frequency extension of the NASA-Ames unsteady, small- 
disturbance code LTRAN2, provides more accurate unsteady results as evi- 
denced by experimental comparisons. The modified code is a versatile tool 
capable of performing reasonably accurate inviscid calculations in both 
linear and nonlinear flow regimes. Results from the improved code may be 
obtained at no extra expense to the user. LTRAN2-HI has now become the 
default option of the NASA-Ames code LTRAN2. 
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□ LTRAN2 (LINEAR) 

A LTRAN2-HI (LINEAR) 

LINEAR THEORY 




Cg o;^ ICg I sin (cot ~ 0g) 


(a) Real and imaginary components. (b) Amplitude and phase. 

Iigure 1.- Lift coefficients vs reduced frequency for a pitching flat plate, 

Mo, 0.7, k = wc/U^, ct - olq ct^ sin tot. 
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Cfi - [Rej) sin cjt - lni|^ cos wtj 

( 3 ) Re3 1 and imaginary 
Figure 4.- Lift 


□ LTRAN2 
A LTRAN2-HI 

© experiment 
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Cg = ICjj^lsin (tut - (p^) 
phase. 

a pitching NACA 64A010 
I sin cot. 


components. (b) Amplitude and 

coefficients vs reduced frequency for 
airfoil, = 0 . 8 . k = ,oc/U„. a = + , 
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□ LTRAN2 
A LTRAN2-HI 
© EXPERIMENT 



“ "1 ' 5*n (cot - 0^1 


(a) Ki.*al and inujglnary components. 


(b) Amplitude and pliase. 


Hgure 5.- l.eadlng-edge moment coefficients vs reduced freuuencv for i 
h.« NACA AAAOlO ..IrfoU. H,. . 0.8. . . ,„c/U„, T.'Mr,;, 
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Hgure 6.- Unsteady lower-surface pressure coefficients ( 
components) for^a pitching NACA 64A010 airfoil, - 0 



Cp ■ (Re sin tjt - Im cos cjt) 




(a) Real and imaginary components. 

Figure 7.- Unsteady lower-surface pressure coefficients for a 
NACA 6AA010 airfoil, - O.U. k n ,.,c/U^ - 0.1. 
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Cp * [Re $ln a»t - Im cos wt] 




Figure 8.- Unsteady lower-surface pressure coefficients (real and imaginary 
components) for a pitching NACA 64A010 airfoil, M„«0.8, k = uc/U„ = 0.2. 
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Figure 9.- Unsteady lower-snrfa^.. 
co»po„e„,., ror puching 


and imaginary 
<^c/Uoo = 0.3. 
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Cp » a-j (Re sin cjt - Im cos cjt] 




Figure 10.- Unsteady lower-surface pressure coefficients (real and imaginary 
components) for a pitching NACA 64A010 airfoil, ® 0.8, k. = (.ic/Uoo = 0.4. 
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